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Abstract. 

We predict a dynamical classical superfluid-insulator transition in a Bose-Einstein 
condensate (BEC) trapped in a combined optical and axially-symmetric harmonic 
potentials initiated by a periodic modulation of the radial trapping potential. The 
transition is marked by a loss of phase coherence in the BEC and a subsequent 
destruction of the interference pattern upon free expansion. For a weak modulation of 
the radial potential the phase coherence is maintained. For a stronger modulation and a 
longer time of holding in the modulated trap, the phase coherence is destroyed signaling 
a classical superfluid-insulator transition. The results are illustrated by a complete 
numerical solution of the axially-symmetric mean-field Gross-Pitaevskii equation for a 
repulsive BEC. Suggestion for future experiment is made. 
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The experimental loading of a cigar-shaped Bose-Einstein condensate (BEC) in 
both one- [T||2^ and three-dimensional [3j optical lattice potentials has allowed the study 
of quantum phase effects on a macroscopic scale such as interference of matter waves [3]. 
There have been several theoretical studies on different aspects of a BEC in one- jS] 
and three-dimensional jO] optical lattice potentials. The phase coherence between 
different sites of a trapped BEC on an optical lattice has been established in recent 
experiments PJ210III1IE] through the formation of distinct interference pattern when the 
traps are removed. In a one-dimensional optical lattice potential the expanding pattern 
consists of a large central piece and two smaller ones moving in opposite directions on 
a straight line In a two-dimensional optical lattice potential the pattern consists 
of a large central piece and eight others on the sides of an expanding square jH]. In a 
three-dimensional optical lattice potential the pattern consists of a large central piece 
and twenty six others on the surface of an expanding cube jSHEj- 

The interference pattern is a consequence of phase coherence in the BEC generated 
due to free quantum tunneling of atoms from one optical lattice site to another 
originating in the superffuid state of the system [HUH]. Equal phase at all points or 
a slowly (and orderly) varying phase are the ideal examples of coherent phase. On the 
other hand, a rapidly (or arbitrarily) varying phase in space is usually incoherent. 

It has been demonstrated for a three-dimensional optical trap potential by Greiner 
et al. |3j that, as the strength of the optical potential traps is increased, the quantum 
tunneling of condensed atoms from one optical site to another is stopped resulting in 
a loss of superfluidity and phase coherence in the BEC. Consequently, no interference 
pattern is formed upon free expansion of such a BEC which is termed a Mott insulator 
state. This phenomenon represents a superfluid to Mott insulator quantum phase 
transition. The phase on an optical lattice site and the number of atoms in that site play 
the roles of conjugate variables obeying the Heisenberg uncertainty principle of quantum 
mechanics jO]. In the superfluid state the coherent phase is considered to be known and 
consequently the number of atoms on each site is unknown thus allowing a free movement 
of atoms from one site to another [3]. In the Mott insulator state the phase is entirely 
arbitrary across the optical lattice sites and the number of atoms at each site is fixed 
and their free passage from one site to another is stopped. As the strength of the optical 
potential traps in the Mott insulator state is reduced the superfluidity is restored in a 
short time via a Mott insulator to superfluid quantum phase transition. This reversible 
quantum phase transition may occur at absolute zero (0 K) and is driven by Heisenberg's 
uncertainty principle [3] and not by thermal fluctuations involving energy as in a classical 
phase transition. As the temperature approaches absolute zero all thermal fluctuations 
die out and at K classical phase transitions are necessarily excluded. 

Following a suggestion by Smerzi et al. [TUj, Cataliotti et al. jTTj have demonstrated 
in a novel experiment the loss of phase coherence and superfluidity in a BEC trapped 
in a one-dimensional optical-lattice and harmonic potentials when the center of the 
harmonic potential is suddenly displaced along the optical lattice through a distance 
larger than a critical value. Then a modulational instability takes place in the BEC and 
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it cannot reorganize itself quickly enough and the phase coherence and superfluidity of 
the BEC are lost. The loss of superfluidity is manifested in the destruction of the 
interference pattern upon free expansion. However, for displacements smaller than 
the critical distance the BEC can reorganize itself and the superfluidity is maintained 
[7|lllj. Distinct from the quantum phase transition observed by Greiner et al. [3 j5 this 
modulational instability responsible for the superfluid-insulator transition is classical in 
nature [Tm illj. This process is also different from the Landau dissipation mechanism 
[T01H2] , occurring when the fluid velocity is greater than local speed of sound. When 
Landau instability occurs, the system lowers energy by emitting phonons [T2]. The 
present classical dynamical transition can be well described [r3* |ll()|ll2j by the mean-field 
Gross-Pitaevskii (GP) equation [T3] . 

The above modulational instability is not the unique dynamical classical process 
leading to a superfluid-insulator transition. Many other classical processes leading to a 
rapid movement in the condensate can lead to such a transition The movement 

should be rapid enough so that the BEC cannot reorganize itself to evolve through phase 
coherent states. In [TJ] a rapid translation of the BEC through the optical lattice sites 
leads to the destruction of phase coherence. Here we suggest that a rapid oscillation of 
the BEC may also lead to a superfluid-insulator transition. The oscillation is initiated 
by a periodic modulation of the magnetic trapping potential ~ uo 2 in the radial direction 
via uj 2 — > oj 2 (1 + Asm(Qr)) where r is time, A an amplitude, uj is the radial trapping 
frequency, and Q is the frequency of modulation. Such modulation of the trapping 
potential is known to generate resonant (collective) excitations in the BEC which have 
been studied both theoretically JH] and experimentally in the absence of an optical 
lattice potential. The study of such excitations in the presence of an optical lattice 
potential has just began [18J. Similar collective excitation generated by a periodic 
modulation of the atomic scattering length [TJ] has also been shown to lead to a classical 
superfluid-insulator transition [20J. 

In the quantum phase transition |3] the Mott insulator state has a perfectly smooth 
probability distribution (modulus of the wave function) across the optical lattice sites 
whereas the phase of the wave function across the optical lattice sites remains entirely 
arbitrary. In the dynamical classical transition considered in this work, because of 
classical oscillation of the BEC, the insulator state in the joint traps is marked by a 
partially disturbed (nonsmooth) probability distribution across the optical lattice sites in 
addition to the loss of phase coherence. However, the information about the destruction 
of superfluidity in both quantum and classical cases is not solely contained in the initial 
probability distribution. Consequently, the BEC needs to be released from the joint 
traps and the formation of the interference pattern studied for a definite conclusion 
about the destruction of superfluidity. 
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As the present transition is classical or mean-field-type in nature, we base the 
present study on the numerical solution of the time-dependent mean-field axially- 
symmetric GP equation [Q] in the presence of a combined harmonic and optical 
potential traps. The time-dependent BEC wave function \&(r;t) at position r and time 
t is described by the following mean-field nonlinear GP equation [T4] 



— ih 



d h 2 V 2 



+ V(r) + gN\^(r;t)\' 



tf(r;t) = 0, 



dt 2m 

where m is the mass and N the number of atoms in the condensate, g 



(1) 

ATrh 2 a/m the 

strength of interatomic interaction, with a the atomic scattering length. In the presence 
of the combined axially-symmetric and optical lattice traps V(r) = \mu 2 (p 2 + i> 2 z 2 ) + 
Vqp^ where uj is the angular frequency of the harmonic trap in the radial direction p, 
vuj that in the axial direction z, with v the aspect ratio, and V Q p^ is the optical lattice 
trap introduced later. The normalization condition is f dr\ty(r; t)\ 2 — 1. 

In the axially-symmetric configuration, the wave function can be written as 
^(r, t) = ip(p, z, t). Now transforming to dimensionless variables p = \/2p/l, z = \/2z/l, 

r = tu>, I = y/h/ (mcu), and (p(p, z; r) = p\J I 3 / \/8ip(p, z;t), (JT]) becomes 
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tp(p,z;r) = 0, 
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Na/l. In terms of the one-dimensional probability P(z,t) = 2ir 



f °° dp\ip(p, z, r)\ 2 /p, the normalization of the wave function is given by dzP(z, t) = 
1. The probability P(z, t) is useful in the study of the present problem under the action of 
the optical lattice potential, specially in the investigation of the formation and evolution 
of the interference pattern after the removal of the trapping potentials. 

In the experiment of Cataliotti et al. [7] with repulsive 87 Rb atoms in the hyperfine 
state F = l,mp = —1, the radial trap frequency was u = 2tt x 92 Hz. The 
optical potential created with the standing-wave laser field of wavelength A = 795 
nm is given by = VoE r cos 2 (kiz), with Er = h 2 k 2 L /{2m) ) ki = 2ti/\, and Vo 

(< 12) the strength. For the mass m = 1.441 x 10~ 25 kg of 87 Rb the harmonic 
oscillator length I = y/KJ (moj) = 1.126 pm and and the present dimensionless time 
unit u^ 1 = 1 / (27r x 92) s = 1.73 ms. In terms of the dimensionless laser wave length 
A = \/2X/l ~ 1, the dimensionless standing-wave energy parameter E R /(ftw) = 47r 2 /Ag. 



Hence in dimensionless unit V^p^ of (j2J is 



opt 



Va 



47T 2 



cos 



2tt 



(3) 



Although we employ the dimensionless space units p and z and time unit r in numerical 
calculation, the results are reported in actual units r /im, z pm and t ms. In the 
conversion we used the parameters of the experiment of Cataliotti et al. [7j, e.g., p = 0.8p 
pm, z = 0.8z pm, and t = 1.73r ms. 



Super fluid- Insulator Transition in a Condensate 



5 



We solve (J2J) numerically using a split-step time-iteration method with the Crank- 
Nicholson discretization scheme described recently [221 • The time iteration is started 
with the harmonic oscillator solution of (J2J) with n = 0: (p(p,z) = [^/(87T 3 )] ly/4 
p e -(p 2 +^ 2 )/4 j2J] The nonlinearity n and the optical lattice potential parameter Vq are 
slowly increased by equal amounts in lOOOOn steps of time iteration until the desired 
value of n and Vq are attained. Then, without changing any parameter, the solution so 
obtained is iterated 50 000 times so that a stable solution is obtained independent of 
the initial input and time and space steps. 

The one-dimensional pattern of BEC on the optical lattice for a specific nonlinearity 
and the interference pattern upon free expansion of such a BEC have been recently 
studied using the numerical solution of (J2J) [13J. Here we study the destruction of this 
interference pattern after the application of a periodic modulation of the radial trapping 
potential in (j2J) via 

p 2 /4^(/3 2 /4)[l + Asin(fir)], (4) 

while the axial trapping potential is left unchanged [TT)]ll7j. In the present model study 
we employ nonlinearity n — 5, the axial trap parameter v = 0.5, and the optical lattice 
strength Vo = 6 throughout. First we calculate the ground-state wave function in the 
combined harmonic and optical lattice potentials. 

Modulation (0J) of the trapping potential may lead to resonant oscillation of the 
BEC and such resonances have been studied in the case of harmonic trap alone [IB]. 
For a very small A, prominent resonances appear in the BEC oscillation when the 
modulation frequency Q is an integral multiple TV of the harmonic oscillator frequency 
uj (O = MuS) [in]- This is quite expected from our wisdom in linear classical physics 
where resonances appear when the driving frequency is a multiple of the characteristic 
frequency of oscillation. Actually, for a finite A, resonances appear for a band of 
modulation frequency [XBj = Muj ± Ajy) where Ajy defines the spread of frequency 
values. The resonance becomes more prominent with the increase of the parameter 
A or A/". At resonance the BEC executes rapid oscillation which is responsible for the 
destruction of superfluidity in the BEC via a classical dynamical transition provided that 
the time of stay of the BEC in the modulated magnetic trap, called hold time, is larger 
than a critical value. A large value of A and/or Q facilitates the dynamical classical 
superfluid-insulator transition. We illustrate this fact in the following for Af = 1 and 2. 
A careful study of the present phenomenon will aid in the understanding of resonance 
in nonlinear physics. Although the generation of resonance in different linear problems 
is well understood, the same in nonlinear physics is just starting. 

As the present calculation is performed with the full wave function without 
approximation, phase coherence among different wells of the optical lattice is 
automatically guaranteed in the initial state. As a result when the condensate is 
released from the combined trap, a matter- wave interference pattern is formed in a few 
milliseconds as described in [T3] • The atom cloud released from one lattice site expand, 
and overlap and interfere with atom clouds from neighboring sites to form the robust 
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Figure 1. Onc-dimcnsional probability P(z,t) vs. z and t for the BEC on optical 
lattice under the action of modulation J3J with = uj and A = 0.5 in the radial 
magnetic trap and upon the removal of the combined traps after hold times (a) 35 ms, 
(b) 52 ms, and (c) 69 ms, in the modulated radial trap. 

interference pattern. The pattern consists of a central peak and two symmetrically 
spaced peaks, each containing about 10% of total number of atoms moving in opposite 
directions jjj. 

Next we consider an oscillating BEC in the combined harmonic and optical traps. 
If we introduce the modulation of the radial trapping potential (jlj) after the formation 
of the BEC in the combined trap, the condensate will be out of equilibrium and start 
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to oscillate. As the height of the potential-well barriers on the optical lattice is much 
larger than the energy of the system, the atoms in the condensate will move by tunneling 
through the potential barriers. This fluctuating transfer of Rb atoms across the potential 
barriers is due to Josephson effect in a neutral quantum liquid J7| . We demonstrate that 
the phase coherence between different wells of the condensate can be destroyed during 
this process for rapid oscillations with large amplitude and/or frequency and no matter- 
wave interference pattern will be formed after the removal of the joint traps. 




Figure 2. One-dimensional probability P(z,t) vs. z and t for the BEC on optical 
lattice under the action of modulation Q with il = 2u and A — 0.5 in the radial 
magnetic trap and upon the removal of the combined traps after hold times (a) 17 ms, 
and (b) 35 ms in the modulated radial trap. 

Now we explicitly study the destruction of superfluidity in the condensate upon the 
application of modulation when the BEC is allowed to stay in this modulated trap for 
a certain interval of time more than a critical value (hold time). For small A and Q away 
from resonance, the BEC executes slow oscillation maintaining the phase coherence. For 
large A and Q and near resonance, the BEC executes rapid oscillation [T7|IT8] which 
results in a destruction of superfluidity. The destruction of superfluidity for a larger hold 
time in the modulated trap manifests in the disappearance of the interference pattern 
upon free expansion which can be studied experimentally. 

For numerical simulation we allow the BEC to evolve on a lattice with r < 15 
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/j,m and 20 /zm > z > —20 /im after the modulation (p|J) with f2 = u; and A = 0.5 is 
applied and study the the system after different hold times. The probability densities 
P{z,t) are plotted in figures 1 (a), (b) and (c), for hold times 35 ms, 52 ms, and 69 ms, 
respectively. For the hold time of 35 ms prominent interference pattern is formed upon 
free expansion as we can see in figure 1 (a). The interference pattern is slowly destroyed 
as hold time is increased as we can see in figures 1 (b) and (c). In figure 1(a) three 
separate pieces in the interference pattern corresponding to three distinct trails can be 
identified. However, as the hold time in the displaced trap increases the maxima of the 
interference pattern mixes up and finally for the hold time of 69 ms the interference 
pattern is completely destroyed as we find in figure 1 (c). 

As the BEC is allowed to evolve for a substantial interval of time after the 
modulation (JH) of the radial trapping potential is applied, a dynamical instability 
of classical nature sets in which destroys superfluidity jHH [TT] . This has been 
explicitly demonstrated in the present simulation which results in the destruction of 
the interference pattern. 

The destruction of superfluidity is facilitated as the amplitude A or frequency Q 
of the radial modulation is increased. We demonstrate this for an increase in Q in 
the following. In figures 2 (a) and (b) we present the evolution of probability P(z, t) 
after the application of the modulation with A = 0.5 and Q = 2uo. With the increase 
of Q from uj to 2u, the destruction of superfluidity is facilitated as one can find from 
figures 2. The superfluidity is destroyed for a hold time of 35 ms for Q = 2uj (figure 
2), whereas it is maintained for the same hold time for Q = u (figure 1). An increase 
in the value of the parameter A also increases the resonant oscillation and we verified 
that the destruction of superfluidity is also facilitated in the process. However, we do 
not present that study here. 

In conclusion, using the explicit numerical solution of the GP equation we have 
studied in detail the destruction of superfluidity in a cigar-shaped condensate loaded in 
a combined axially-symmetric harmonic and optical lattice traps upon the application 
of a modulation of the radial trapping potential near resonance. In the absence of 
modulation, the formation of the interference pattern upon the removal of the combined 
traps clearly demonstrates the phase coherence jHl^j. The superfluidity is maintained 
for a slow modulation (J3J) of the radial trapping potential away from resonance. For rapid 
modulation of large amplitude and/or frequency near resonance, there is a superfluid- 
insulator classical dynamical transition, provided that the BEC is kept in the modulated 
trap for a certain hold time. Consequently, after release from the combined trap 
no interference pattern is formed. For smaller amplitude A and/or frequency Q of 
modulation (J1J near resonance, the superfluid-insulator transition occurs for a larger 
hold time and vice versa. It is possible to study this novel phenomenon experimentally 
and compare with the present theoretical prediction. 
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